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ABSTRACT 

Strong magnetic fields play a crucial role in the removal of angular momentum from collaps- 
ing clouds and protostellar discs and are necessary for the formation of disc winds as well as 
jets from the inner disc and indeed, strong large-scale poloidal magnetic fields are observed 
in protostellar discs at all radii down to ~ IORq. Nevertheless, by the time the star is visible 
virtually all of the original magnetic flux has vanished. I explore mechanisms for removing 
this flux during the formation of the protostar once it is magnetically disconnected from the 
parent cloud, looking at both radiative and convective protostars. This includes a numerical 
investigation of buoyant magnetic field removal from convective stars. It is found that if the 
star goes through a fully convective phase all remaining flux can easily be removed from the 
protostar, essentially on an Alfven timescale. If on the other hand the protostar has no fully 
convective phase then some flux can be retained, the quantity depending on the net magnetic 
helicity, which is probably quite small. Only some fraction of this flux is visible at the stellar 
surface. I also look at how the same mechanisms could prevent flux from accreting onto the 
star at all, meaning that mass would only accrete as fast as it is able to slip past the flux. 
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1 INTRODUCTION 

Stars form from clouds in which gravitational, thermal, rotational 
and magnetic energies are generally comparable to one another. 
Once a star has formed, the thermal and gravitational energies are 
still comparable to each other but the magnetic energy (as inferred 
from Zeeman measurements) is at least six orders of magnitude 
lower in even the most strongly magnetised stars. 

This is not what one prima facie expects if a cloud retains 
its magnetic flux as it collapses, since gravitational and magnetic 
energy have the same scaling with radius, i.e. £ a 1 / R. In con- 
trast, thermal and rotational energy rise faster, and must therefore 
be continually extracted from the cloud to allow its collapse. Ex- 
traction of thermal energy happens via radiation and represents no 
barrier to star formation in theory or otherwise (at least for stars 
under about lOM©); extraction of rotational energy is less well un- 
derstood but theoretical mechanisms exist, using magnetic fields or 
gravitational i nstability, and there is plenty of evidence for mag- 
netic braking {Gillis et alii 19741 . 1 19791 : iMouschovias & Paleo logou 
Il979l : lstahler & Pallall2005l and refs. therein). 

The relative strength of the gravitational and magnetic fields 
is often expressed as a dimensionless mass-to-flux ratio, defined as 
A = M/$, where M and <J> are the mass and magnetic 

flux, or locally in a disc context as 27rG 1//2 E/ B z where B z and 
E are the field normal to the disc, and surface density; this ratio 
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is conserved if flux freezing is valid. It is related to the gravita- 
tional and magnetic energies by (ignoring factors of order unity) 
A 2 = |-E g rav|/-Emag- A cloud with A > 1 is said to be 'magneti- 
cally supercritical' and will collapse, in the absence of significant 
thermal or rotational energy. Conversely a cloud with A < 1 is 
'magnetically subcritical' and the magnetic field supports the cloud 
against gravity. 



Observations of cloud cores and the ISM (e .g. Crutcher et al.l 
1 19991 : iHeiles & Trolandll2005l : iGirart et alj|2006l . l2009h show that 
clouds do contain strong magnetic fields and appear to be mildly 
magnetically subcritical on scales above ~ 1000 AU. The rota- 
tional and thermal energies appear to be relatively small on these 
large scales. The so-called "magnetic flux problem" can be di- 
vided into two parts, the first of which is: how is sufficient mag- 
netic flux lost so that the cloud becomes magnetically supercriti- 
cal? Gravitational cont raction down to that sc ale could proceed via 
ambipolar diffusion i Mestel & Spitzerll 19560 or some other diffu- 
sive process, and dynamical collapse begins once supercriticality 
is reached. Recent obs ervational support for this model is given by 
iDavidsoneTall l bOllh . 

The second part of the magnetic flux problem - that ad- 
dressed here - is how to remove the remaining flux once the cloud 
has become supercriti cal, to explain the incredibly weak magnetic 
fields seen in stars (e.g. lMestel & Spitzerll956l : lNakanoll984l : iGallil 
2009). The supercritical collapse might take place too quickly to 
allow ambipolar diffusion to remove the rest of the flux before 
flux freezing resumes owing to rising ionisation fraction (although 
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opinio n s differ on this point, see 10 dl998h ; l5esch & Mouschoviasl 
(2001); Ban eriee & Pudrit j d2006h amongst others). 

In this paper I look at a solution to this problem, namely that 
the collapse can proceed with as little or as much loss of flux as dis- 
sipative processes allow, and that the excess flux is destroyed once 
the protostar forms, or rather, once it has become disconnected in 
some sense from its parent cloud. In the next section I review rel- 
evant observational and theoretical results of collapse from 1000 
AU to protostar formation. In section [3] I explore the mechanisms 
for destroying flux in stars via magnetohydrodynamic (MHD) in- 
stability in non-convective stars and buoyancy in convective stars, 
before investigating the latter numerically in section U.2.21 1 discuss 
the results in section|4]and conclude in section|5] 



2 ACCRETION AND PROTOSTAR FORMATION 

I now summarise relevant results and evidence that at least part 
of the solution to the magnetic flux problem must lie in, or in the 
immediate vicinity of, the protostar. 



2.1 Inwards advection of magnetic flux in a disc 

Once a cloud has become supercritical, it can collapse dynamically. 
Magnetic braking becomes ineffective once the collapse is super- 
Alfvenic, so the rotational energy becomes larger in relation to the 
other energies. Normally this leads to formation of a disc of ra- 
dius 100 — 1000 AU, although some system s might lack a disc 
dStassun et all 19991 1200 ll iRebull et alj|2006f) . However, we shall 
look here at the magnetic flux accreted via a disc. 

Gas in an accretion disc can lose its angular momentum 
and spi ral inwards by passing it either vertically to a disc wind 
or jet jBlandford & Pavna Il982h . or radially outwards to disc 
gas exterior to itself via some instabil ity-driven turbulence. The 
instability could be gravitational (e.g. Ijappsen & Klessenl 120041 



Vorobvov & Basul2006l:lzhu et alj2009h . magnetorotational (MRI 
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-. .Pessahl l2010l and refs. therein) perhaps 
with buoyancy (e.g. iKeppens et al. [ |2002l) or purely hydrodynam- 
ical dLithwickl2009l) . 

Jets are observed in many objects, are launc hed from radii ~ 1 
AU ( iBacciotti et alj|2002l ; lAnderson et alj|2003f) and often appear 
together with a slower, le ss collimated wind. Around a tenth of th e 
accreted mass is ejected dKurosawa et al. 2006; Podio et al]|2006t) . 
There is evidence in favour of both theories; for instance the out- 
flow model is supported by measurements of the angular momen- 
tum of jets dRav et al J 2007). but not by the fact that discs are self- 
luminous, since an outflow which carries all of the angular momen- 
tum also removes all of the accretion energy. It is likely therefore 
that discs are outflow- or turbulence-dominate d in different zones, 
at different times and in different objects (see lCombet et alj|2010l 
and refs. therein for a recent comparison of the two). 

There is strong evidence that discs do contain str ong, ordered, 
net poloidal flux. This is crucial for the o utflow mo del ( Ouv ed et al.l 
1 19971; iBaneriee & Pudrit j l200d. 120071 ; iBeckwith et alj|2008l) but 
possibly also helpful in increasing the efficiency of the MRI, bring- 
ing the Shakura-Sunyaev a paramet er produced in simulations up 
to the observationally-inferred value dKing e t al. 2007; Pessah et al.l 
l2007h . Moreover, there are direct measureme nts of the mag- 
netic fi eld in protostellar discs at various radii. Iviemmings et al.l 
d2010h used methanol masers in the massive-protostar system 
CepheusAHW2, finding a field of 23 mG at a disc radius of ~ 1000 



AU and a local mass-to-flux ratio of A ~ 1.7. The field is large- 
scale and has a di rection similar t o that in the surrounding cloud. 
Using meteorites, iLevv & Sonett ( 1978) measured 3 G at radius 
1 AU in the pro to-solar system, tr anslating to A ~ 5, using ob- 
servations (e.g IWilner & Lay||2000h giving surface densities of or- 
der 10 4 g cm" 2 at this radius. Finally, in the FU Orionis system 



iDonati et all d2005l) used the Zeeman effect to measure a vertical 
field component of 1 kG with a filling factor of 20% at a radius of 
0.05 AU, i.e. KLR©. With a surface density of ~ 4 x 10 4 g cm" 2 , 
this gives a mean mass-to-flux ratio A ~ O.10This poloidal field 
is accompanied by a somewhat weaker toroidal component, con- 
trasting with local dynamo models where a predominantly toroidal 
field is generated in situ. Although these observational results are 
approximate and relate to different radii in different objects, we see 
that discs are close to magnetically critical even very close to the 
protostar - there is no evidence for A anywhere near as high as in 
stars (A = 10 3 - 10 s , see below). 

It is important to note that a small A(r) = 
27rG 1//2 £(r)/ B z (r) in the disc does not necessarily mean 
that mass and flux are being adverted inwards in that same 
ratio (and if the aforementioned value of A ~ 0.1 is true, this 
is impossible). The mass might be slipping inwards past the 
magnetic field lines; all we can say with certainty is that the flux 
present in the disc was adverted inwards at some point in the 
past. In the turbulent-disc scenario, Ivan Ballegooijen ( 1989) used 
geometrical arguments to show that the turbulent diffusion should 
cause the magnetic field to diffuse outwards relative to the gas at 
the roughly same speed at which the gas moves inwards . There 
may be ways around this: for instance, Spruit & Uzde nskv ( 2005) 
proposed a mechanism to advert flux inwards in discrete clumps, 
which is perhaps s upported by the aforementioned observations of 
IDonati et al.l d2005h . In contrast, in the disc-wind/jet model, one 
expects prima facie that flux is adverted inwards since accretion is 
fast and there is no turbulent diffusion. 

If in the steady state the star is accreting mass and flux in the 
ratio A,, then mass and flux must be passing through each sur- 
face of constant radius in the disc in the same ratio (ignoring out- 
flow), even though the local ratio will in general be much lower 
A(r) <C A», requiring almost perfect slippage at all radii. Since 
the properties of the disc vary significantly over the large range in 
radii, it is natural to infer that there is some feedback mechanism 
which prevents flux from moving inwards faster than it can be ab- 
sorbed into the star, which in turn will limit the mass accretion rate 
to that at which the mass can effectively slip past the field lines. 
This mechanism could work locally, 'feeling' some quantity such 
as the radial gradient in field strength, or it could be cyclic. Since 
observations show that A(1OJ?0) <C A» it seems likely that this 
feedback originates from a bottleneck in the central region where 
it is noteworthy that a high ionisation fraction renders ambipolar 
and Ohmic diffusion ineffective. Finally, note that we do not know 
the value of A, while an embedded star is accreting - it may well 
be very much lower than the mass-to-flux ratio of stars which have 
finished accreting. It is possible that much of the flux is lost once 
the main accretion phase is over and the star has become detached 
from its surroundings. 

The purpose of this paper is to show that any excess flux can 
be destroyed once the protostar forms. Excess flux is taken to mean 
the difference between the flux accreted onto the protostar and the 



1 Note that j3 ~ A 2 (/i/r)M*/M disc where /3 = 8ttP/B 2 ; in a disc 
therefore it is possible to have both /3 > 1 and A < 1. 
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flux later observed on the star. The easiest way to imagine how this 
can happen is to assume first that all flux which survives diffusive 
processes during the supercritical core collapse is accreted onto the 
protostar, which therefore at least initially contains the entire ac- 
creted flux. The star then becomes magnetically 'detached' from 
its parent cloud, and then after accretion has finished the flux is de- 
stroyed according to the mechanisms described in section [3] One 
might actually expect these mechanisms to work continuously as 
accretion is ongoing, so that the flux of the infalling material is lost 
as soon as it arrives at the protostar; this scenario is discussed in 
section [33l 

2.2 Conditions in the proto- and main-sequence star 

Of crucial importance for the following discussion is the appear- 
ance of convection, which appears in protostars when a tempera- 
ture of 10 6 K is reached and deuterium fusion begins. If the mass 
remains below OAMq, the star remains on the fully-convective 
Hayashi track all the way onto the main-sequence (MS) when hy- 
drogen ignites. Above OAMq a radiative interior eventually de- 
velops, either after accretion has ceased, whereupon the star e nters 
the Henyey track dHenvev et ai1ll955l ; IPalla & Stahlej|l993h and 
moves to the left on the HR diagram, or in a more massive protostar 
while accretion is still ongoing - deuterium burns in the convective 
shell at the same rate at which it is accreted. Ab ove about 4Mp) 
conve ction retreats right to the surface of the star jPalla & StahleJ 
1 1993b and may dissappear entirely or have little effect on the ac- 
creting gas. This marks a fundamental difference between high- 
mass and low-mass stars, with the all the material in stars below 
AMq having experienced convection, but only the inner 4Mq of 
more massive stars having experienced convection. In addition, it 
is possible that very massive protostars somehow bypass the fully 
convective phase during growth up to 2Mq, retaining a radiative 
interior at all times. 

On the MS, stars below OAMq, which are fully convective, 
and stars between 0.4 and 1.5Mq, which have a convective enve- 
lope, display fluctuating dynamo fields. It is not obvious from first 
principles whether a magnetic-convective steady state should de- 
pend on the initial conditions; all we can say is that no such depen- 
dence is apparent from the observations, which do however show 
a clear positive correlation b etween r otation speed and magnetic 
activity jPizzolato et al ]2003l ; see also lMorin et al]20ld) . 

In contrast, in MS stars with a radiative envelope (> 1.5M©) 
the observations of magnetic fields are the following. Of t he less 
massi ve stars (up to 6Mq) some fraction (~ 10%, see e.g. lPowei] 
120071) display a weak, large-scale field of 0.2 — 30 kG; the so- 
called Ap and Bp stars. At higher masse s (O and early B) it 
seems that a similar fraction is magnetic jGnmhutetal7l201ll) . 
The fields strengths measured correspond to a mass-to-flux ratio 
of A < 10~ 3 assuming similar interior and surface field strengths. 
These fields are presumably in equilibrium, since there is no possi- 
bility of a contemporary dynamo. The rest of the p opulation, how- 
ever, s eem to have field s at the gauss-level or lower (Lignier es et al.l 
l2009l ; |Petitetai]|201ll) , which may be for instance produced by 
a subsurface-convective-dynamo mechanism (Cantiello & Braith- 
waite, in prep.) or be evolving decaying dynamically on a timescale 
given by t a /P where ta and P are the Alfven timescale and ro- 
tation period (Braithwaite & Cantiello, in prep.) In any case, the 
important point is that in contrast to convective stars, the mag- 
netic properties of radiative stars do depend on the initial condi- 
tio ns. For a recent review o f magnetism in main-sequence stars 
see lDonati & Lan dstreet 2009. Fields can also be observed directly 




Figure 1. The magnetic disconnection of a collapsing core/protostar from 
the rest of the cloud (for clarity, very simplistic and not to scale), (a) A col- 
lapsing overdense region draws field lines together, the classical 'hourglass' 
picture, (b) Further collapse, (c) Reconnection occurs outside the protostar, 
which becomes disconnected from the cloud, (d) A super-Alfvenic wind 
ensures that the star is completely causally disconnected from its surround- 
ings, and the flux in the star can now be destroyed. The wind expels the 
surrounding cloud material. 

on stars which are still accreting; in both convective and non- 
convective stars the fields may not be very diffe rent from those on 
the main-sequence (see e.g. lAlecian et al", 2009). 

This agrees with what one expects theoretically. In a convec- 
tive star, the magnetic field reaches a steady state independent of 
the initial conditions. In a radiative star the magnetic field evolves 
on a dynamic timescale until it finds some stable equilibrium, the 
geometry and strength of which does depend on the initial condi- 
tions. 



3 DESTRUCTION OF FLUX IN A STAR 

Once accretion has stopped and the star has developed a super- 
Alfvenic wind, the star is 'detached' from the parent cloud and 
its flux no longer needs to be conserved. In fact, this detachment 
may happen earlier than the complete cessation of accretion, when 
the accretion rate decreases, and the star could possibly detach and 
reattach many times before final detachment, in response to FU On- 
onis accretion episodes. A possible sequence of events is illustrated 
in fig.Q] and although the actual path between stages (a) and (d) is 
uncertain, the detail is unimportant. I now look at two mechanisms 
to destroy flux in the idealised case of a fully detached star, which 
work respectively via MHD instability and reconnection to equi- 
librium inside a radiative star (section [3~Tt and buoyant expulsion 
from the interior of a convective star (section [3^21 . Afterwards, in 
section [3~3l I look at the more realistic scenario where these mech- 
anisms proceed while the star is still accreting. 

3.1 Non-convective stars 

In the simplest picture, a stably-stratified star initially contains a 
symmetric, dipole-like purely-poloidal field which was inherited 
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conte xts from the laborat ory dChui & M offatt 199H; lHsu & Bellanl 
l2002t) to the solar corona fchang & Lowl2003l) . Therefore: 



Figure 2. Form of the instability in a poloidal field. On the left is shown the 
equilibrium field-line loops; on the right the form of the unstable motions 
at some particular azimuthal wavenumber. 



from the cl oud. Any purely p oloidal field i s subject to an MHD 
instability dMarkev & Tavletl 1 19731 1 19741; [Flowers & Rudermanl 
ll977l:lBraithwaitdl2007l : lMarchant et al.ll2010h which is similar to 
the instability in a pair of initially aligned bar magnets which rotate 
until they lie with opposite poles next to each other. The free en- 
ergy is the volume integral of B 2 /8n outside the magnets, which 
drops as a result of the realignment. A fluid star can be thought of 
as a collection of bar magnets; the magnets do not stay aligned and 
the instability grows on the dynamical (Alfven) timescale, around 
10 years in a solar-type star with a field strength of 1 kG - in any 
case much shorter than any formation timescale. See fig [2] Even- 
tually length scales become small enough so that magnetic flux is 
destroyed inside the star, resulting in the complete destruction of 
the magnetic field. Note that even if the field is not purely poloidal 
but is 'wound-up' to some degree, there is nothing to stop it un- 
winding first and then developing this instability. 

If the field is not symmetrical as above but has some net 'twist' 
then MHD instability will not destroy it entirely. Instead, it will 
evolve on the Alfven timescale into some stable equilibrium. The 
main factors which determine the nature of the resulting equilib- 
rium are a) the net twist, or more formally the magnetic helicity 
and b) the initial radial distribution of magnetic flux. 



3.1.1 Helicity 

Let us imagine the relaxation to magnetohydrodynamic equilibrium 
in a star of radius R, which contains a magnetic field of energy 
E = (47r/3)i? 3 B 2 /87r, where B is the r.m.s. magnetic field. In 
the following, quantities initially and at (finally) equilibrium are 
marked with the subscripts i and f respectively. We can say the 
following about the equilibrium state. 

The reconnection destroys magnetic energy on small length 
scales but has little effect on the magnetic helicity, a global quan- 
tity defined as the volume integral of the scalar product of the mag- 
netic field with its vector potential H = (l/87r) JA ■ B dV. It 
can be sho wn that in the case of infinite conductivity, helicity is 
conserved (Wolt jei]|l958l) . Helicity has units of energy x length 
and so is present more in the larger structures than is the energy - 
and it is approximately conserved during reconnection taking place 
on small scales, a property which has been very useful in many 



H f « Hi. 

Consideration of dimensions gives us 

\H f \ =iptR*E t) 



(1) 



(2) 



where ipf is a dimensionless parameter of the equilibrium; we can 
solve for Ef and Ht once we know its value. In a large-scale equi- 
librium ipi should be of order unity; in a small-scale equilibrium it 
will be smaller but note that during a relaxation to equilibrium at 
constant helicity an equilibrium with higher ip{ is always favoured 
since it represents a lower energy state. 

So, the energy and strength of the magnetic field on the main 
sequence is determined by the helicity of the parent cloud, which 
can be expressed as |i? c t| = ipciRciEd where ip c \ is the degree of 
twist in the cloud, which can have any value from to order plus or 
minus unity (helicity is either positive or negative depending on the 
sense of the twist). This logic was rece ntly confirmed num erically 
in the context of intergalactic bubbles teraithwaitell2010l) . where 
a relaxation to equilibrium was found approximately to conserve 
helicity, upon which the final equilibrium energy therefore depends. 

Finally, note that in the standard hourglass model, helicity is 
zero and so one expects the field to disappear completely after 
disconnection from the cloud. A non-zero helicity requires some 
asymmetry between the hemispheres. 



3.1.2 Stratification 

Assuming that the magnetic energy is much less than the gravi- 
tational, a positive radial entropy gradient prevents the gas from 
moving significantly in the radial direction (we speak of a stably- 
stratified star), so that during the relaxation to equilibrium the gas 
motion is confined to spherical shells0This means that the absolute 
flux through a spherical shell cannot increase during relaxation to 
equilibrium. This can be seen by considering two or more regions 
on a spherical shell of positive and negative radial component of 
magnetic field B r ; see fig. [3] As the fluid moves around on the 
spherical shell discontinuities (current sheets) develop between re- 
gions of different B r (not just between positive and negative B r 
but between any differing values) and the result of these sheets is 
that gas from either side with initially different B r ends up with 
some B r which is an average of the two. Indeed in general, dis- 
sipative processes can only lead to a drop in the spherical surface 
integral § |B r |dS at any given radius. Therefore, if the star ini- 
tially contains a strong magnetic field in the centre and a weak field 
at the surface, which is what one might expect if each fluid ele- 
ment conserves its flux during accretion and compression so that 



B 



oc p 



2/3 



then the resulting equilibrium will also have only a 



weak field at the surface. It is possible for a star to 'hide' a large 
magnetic energy in its interior but show nothing at the surface. In 
fact, if B oc p 2 / 3 thoughout the star then Ap stars, which have mag- 
netic to thermal energy density ratios at the surface of between 1 
and 100 (i.e. the plasma-/3 is 0.01 to 1), have global Exnsig/\E smv | 
ratios of the about same value. Of course, a global ratio of unity 
or greater is impossible since the star would be gravitationally un- 
bound. However, the lower conductivity near the surface means that 
a star which initially has B oc p 2//3 will undergo diffusion which 

2 In addition, the motion is approximately incompressible, so the velocity 
field has just one degree of freedom. 
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current sheet 



spherical shell 




Figure 3. Reconnection on a spherical shell where a discontinuity appears 
in the magnetic field. The grey arrows represent the velocity field; note that 
the motions are essentially two-dimensional. 



leads to the field in a layer below the surface relaxing to a potential 
field, meaning that the field strength at the surface will be similar to 
that at the bottom of the layer, where the plasma-/? is much greater 
than unity; the B <x p 2//3 relation continues to hold deeper in the 
interior. For instance, in a 2Mq ZAMS star the potential-field layer 
should grow to a depth of ~ 0.05J?* after a time 10 7 yr (calculated 
by equating t to L 2 Jr\ where L is depth and r\ the magnetic dif- 
fusivity at that depth); at that depth the density is ~ 10 4 times its 
photospheric value; if a field strength of 3 kG seen on the surface 
is roughly constant down to that depth and then B oc p 2 ^' A below 
that, the global ratio E ml!lg /\E gISLV | ~5x 10~ 4 . Most of the flux 
would be hidden within the star. Alternatively the field strength in 
the interior may depend much less steeply on density; in fact if one 
imagines building a (radiative) star by accretion, adding spherical 
shells from inside to outside (as seems likely from entropy consid- 
erations), it is difficult to see how the fluid elements could possibly 
retain their original flux without a very unlikely and unstable mag- 
netic field configuration. In any case, there is still potential for the 
field strength deep in the interior to be much greater than that seen 
on the surface. 



3.2 Convective stars 

Alternatively, if at some time the star becomes convective then the 
magnetic field will tend to rise to the surface, its energy being dis- 
sipated by reconnection in the atmosphere. While the magnetic en- 
ergy density is greater than the convective, as might be expected 
in a protostar at least initially, the important condition is not the 
convective motion, but the near-isentropic state which this motion 
maintains. Deviations from uniform entropy are also unimportant 
when the field is above equipartition with the convective kinetic 
energy - buoyancy effects are dominated by the magnetic field. 



3.2.1 Buoyant loss of magnetic flux 

A magnetic field in an ise ntropic star has a tendency to rise towards 
the surface (Reiseneggei 2009, and refs. therein). Any magnetised 
region must be in pressure equilibrium with its surroundings, which 
means that its gas pressure must be lower than that in the surround- 
ings. Since the entropy is constant we have p = p(P), so that the 
region must also have lower density. There is no way to stop buoy- 
ant rise to the surface unless reconnection can occur fast enough 
that a global magnetic equilibrium is reached where buoyant re- 
gions are somehow held down by magnetic tension. It is therefore 
informative to estimate the relevant timescales. 

Consider a body of gas sitting in hydrostatic equilibrium in 



a gravitational potential well. The gas contains a non-equilibrium 
magnetic field of r.m.s. magnitude B with a characteristic length 
scale I, which is initially much smaller than the size of the sys- 
tem R (below, the case with initial I ~ R is looked at). Magnetic 
reconnection occurs at a speed 



avA 



(3) 



where the reconnection speed parameter a is of order unity (in the 
solar corona context a value 0.1 is often assumed) and va is the 
Alfven speed. As reconnection proceeds on a timescale of / / (ova), 
in the absence of other effects the length scale I increases until some 
global equilibrium is reached with / of the order the size of the 
system. However, while this reconnection is proceeding, there are 
buoyant motions due to magnetically-induced density variations, 
which also have a length scale /: imagine the buoyant motion of 
a parcel of gas of size I and density p through its surroundings of 
density po. Matching the buoyancy force to the aerodynamic drag, 
we have 



pol 2 v 2 ~ g\po - p\l 3 



(4) 



where v t is the terminal velocity (upwards or downwards depend- 
ing on the sign of po — p) and g is the gravitational acceleration. 
It can easily be shown that the terminal velocity is achieved after 
the region has risen a distance I. From the equation of hydrostatic 
equilibrium we can reorganise {4} to give 



Vt 



1/2 



Ipo - p\ 



I 



Hp J \ Po 



1/2 



Co, 



(5) 



where H p is the pressure scale height and Co is the sound speed 
in the external medium (recall that Cq ~ Po/ Po). We can assume 



so that the motion is very subsoniqj. Pressure balance 

give 

R 2 R 2 

This ignores the anisotropic nature of magnetic pressure; in reality 
neighbouring blobs will push and pull on each other in some com- 
plex way, but we limit ourselves here to looking just at the general 
tendency for buoyant rise. In an isentropic star we have P = Kp J 



Po 



Po 



7- 



Po 



po 



(7) 



in contrast to a star with an extra thermodynamic degree of free- 
dom where buoyant balance p = po can be achieved despite the 
pressure difference: for instance in a star with an ideal-gas EOS 
the temperature can be adjusted. Using {6]l and Q, ((5) becomes 
(dropping factors of order unity and assuming \B — Bq\ ~ B) 



t't 



1/2 



t'A- 



(8) 



This corresponds to the result o f |Parken ( fl975h for the rise of a flux 
tube in a convective zone where the hydrodynamic force associated 
with the convective motion can be ignored. Comparing this to (£/3j> 

3 In the case of X-ray cavities in galaxy clusters, we have p <C po and I ss 
Hp, so that vt ~ Co. From the apparent absence of shocks in these systems 
we infer that the motion is in fact mildly subsonic and cannot therefore 
solve the cooling-flow problem with shock heating. 

4 It is appropriate here to use an 'isotropic' magnetic pressure equal to one 
third of the energy density, as with any relativistic fluid (see e.g. Br aithwaite] 

boich . 
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we see that the condition for reconnection to proceed faster than 
buoyant motion is 



It seems therefore that the relaxation of an arbitrary magnetic field 
with small initial I to equilibrium as seen in simulations of a sta- 
bly stratified star cannot occur in a star with p — p(P) beyond the 
point where the length scale is some fraction of the scale height. If 
on the other hand the initial conditions have a larger length scale I 
than that in @, the buoyant motion is always faster than the recon- 
nection speed and so it is impossible for the magnetic field to reor- 
ganise itself. Regions of strong field will rise to the surface while 
regions of weak field will sink to the centre; at or above the surface 
of the star (in the low-plasma-/3, force-free regime) there could be 
reconnection and loss of magnetic helicity, eventually resulting in 
complete loss of the magnetic field of the star via a continuous pro- 
cess of magnetic buoyant convection. In principle it seems possible 
that in the absence of heat-driven convection, some field might be 
retained and form a global equilibriurr0; however, an equilibrium 
not only requires non-zero helicity as in the case of stable stratifi- 
cation explored in section [3T| but should eventually be destroyed 
by the convective motions. 

However, how can one justify ignoring the convective motion 
in the above analysis? To suppress convection, a field must not only 
be coherent on scales larger than the convective cel ls (~ H v ) bu t 
also at least at equipartition with the thermal energy (Mestel 1970), 
which we can assume is not the case for Hayashi-track stars. How- 
ever, while the energy in the field is greater than the kinetic en- 
ergy of the convection the convection is sub-Alfvenic and therefore 
makes little difference to the magnetic buoyant motions explored 
above. 

In any case, once the magnetic field has decayed to equipar- 
tition with the convection, the convective motion can no longer be 
ignored, the analysis above becoming invalid. The domain of the 
convective dynamo has been entered, and some steady-state is pre- 
sumably reached with the magnetic energy at most comparable to 
the convection energy. Any 'memory' of the original (large) mag- 
netic flux is presumably erased, and the main-sequence magnetic 
field does not depend on the conditions in the parent cloud. 

3.2.2 Numerical model of an isentropic star 

I now describe a numerical model of flux loss from an isen- 
tropic star. The stellar model and computational setup are simi- 
lar to those described in iBraithwaite & Nordlunj|2006l to where 
the reader is referred for a fuller account of the setup of the 
model ; a brief outline is given here. The code used is the STAGGE R 
CODE dNordlund & Galsgaard|[l99lllGudiksen & Nordlundll2005f) . 
a high-order finite-difference Cartesian MHD code which uses a 
'hyper-diffusion' scheme, a system whereby diffusivities are scaled 
with the length scales present so that badly resolved structure near 
the Nyquist spatial frequency is damped whilst preserving well- 
resolved structure on longer length scales. This, and the high-order 
spatial interpolation and derivatives (sixth order) and time-stepping 



(third order) increase efficiency by giving a low effective diffusiv- 
ity at modest resolution (192 3 here). The code includes Ohmic and 
well as thermal and kinetic diffusion. Using Cartesian coordinates 
avoids problems with singularities and simplifies the boundary con- 
ditions: periodic boundaries are used here. 

The simulations model the star as a self-gravitating ball of 
ideal gas (7 = 5/3) of radius R in hydrostatic equilibrium with 
radial density and pressure profiles obeying the polytrope relation 
P oc p 1+1 / n , with the index n set to 3/2 here to give constant 
en tropy, as opposed to the radia tive-star approximation n — 3 
in IBraithwaite & Nordlunoj 2006. Surrounding the star is a hot, 
poorly-conducting atmosphere, which has the effects both of forc- 
ing the magnetic field to relax into a potential (curl-free) state while 
keeping the Alfven speed outside the star numerically convenient, 
since the density is not too low; in other words it is a simple way 
of modelling a vacuum. Since the field strength drops by a large 
factor during the process being modelled, a field-upping routine is 
employed, which artificiall y amplifies the field to keep th e mag- 
netic energy constant (see Brait hwaite & Nordlundl 2006 for de- 
tails), which greatly reduces the computational demand and, per- 
haps more importantly, keeps the Alfven timescale much shorter 
than the diffusive (Ohmic) timescale. At the beginning, the star is 
given an initially random magnetic field containing energy at all 
length scales down to a limit of a few grid-spacings, and the MHD 
equations are integrated in time to follow the evolution of the field. 

In fig. [4] we can see an example of the evolution of a field in 
this numerical model. At the beginning the field is somewhat con- 
centrated into the interior of the star, and rises towards the surface 
on an Alfven timescale. Another example of this, starting with a 
different randomisation of the initial field, is shown in fig. [5] 

Taking a closer look, at some point in time some recognis- 
able structure becomes visible, in the form of flux tubes lying 
quasi-horizontally around the surface of the star. These resem- 
ble the struc ture of non-a xisymmetric equilibria found in previous 
simulations l lBraithwaitell2008l) when a similarly random, but less 
centrally-concentrated initial field (i.e. with significant flux thread- 
ing the surface at t — 0) is left to evolve in a stably-stratified star - 
horizontally-lying flux tubes arranged in some meandering pattern 
below the surface of the stajfl these tubes evolve on a timescale 
given by the Ohmic timescale and/or a large multiple of the thermal 
timescale (whichever is shorter) under the stellar surface, which 
can of course be significant compared to, or longer than, a main- 
sequence lifetime. In contrast, the similar structures found in the 
current study are evolving on a shorter timescale related to the 
Alfven timescale in the stellar interior, due to the buoyant force 
pushing them upwards into the atmosphere of the star, and the 
Ohmic timescale above the stellar surface, as the tubes decay due 
to (relatively low) finite conductivity. It is not immediately obvious 
what this timescale should be, except that it will be rather small 
compared to the evolution timescale of a protostar. In the simula- 
tions, the exterior has a much lower conductivity than the interior, 
which forces the exterior field to relax to a potential (curl-free) con- 
figuration, but it is impossible to make this effect quantitatively re- 
alistic in the numerical model while dynamic-timescale processes 
are ongoing. In all the simulations the decay timescale is a few 
Alfven timescales. In reality the decay above the surface (in the 



5 The range of available equilibria in isentropic stars, or alternatively in 
stars with a barotropic EOS p = p(P), is more restricted than in non- 
barotropic stars since there is only one degree of freedom in adjusting the 
thermodynamic state to balance the Lorentz force, but their existence cannot 
be ruled out. 



B Note that in the simulations it is difficult to define the precise location 
of the stellar surface, as the transition from isentropic and high electrical 
conductivity to isothermal and low electrical conductivity takes place over 
several grid spacings 
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Figure 4. The evolution of an initially turbulent magnetic field in a non-stably stratified star. The three snapshots are at t/r\ = 0, 3.1 and 7.8 where ta is the 
Alfven timescale. After some Alfven times, the field has bubbled up to the surface and shows no signs of reaching an equilbrium. The blue shading represents 
the stellar surface. 




Figure 5. When the field reaches the surface of the star, familar structures 
become visible. Here, field lines are plotted and the surface of the star is 
coloured according to the value of the radial field component (blue & vio- 
let are negative, green & yellow/orange are positive). This run is similar to 
the one shown in fig. [4] but with a different random initialisation; this snap- 
shot is taken after 12.2 Alfven timescales have passed (corresponding to the 
last frame of fig- [4] Flux tubes lie along the surface in apparently random 
patterns. 

low plasma-/? regime) should (helped by convective motions ab- 
sent from these simulations) proceed on a dynamic timescale in 
localised reconnection zones, as it does in the solar corona. 

In fig.[6]we can see the process in a little more detail - in a time 
animation from a third run with a different random initialisation. 
Several runs were done, with different random initialisations and 
with different initial radial profiles of magnetic energy, but in each 
run the result was essentially the same. 

Broadly speaking, these numerical results confirm the predic- 
tion made in section [3T2l This behaviour is distinct from that of 
the same magnetic field in a stably-stratified star, where there is no 
transport of material in the radial direction (Braithwaite & Spruit 



|2004 lBraithwaitell2008l) and consequently no transport of radial 
flux from the interior to the surface. It seems therefore that an MHD 
equilibrium cannot be reached from these initial conditions in an 
isentropic star; however it is impossible to rule out at this stage that 
different initial conditions, perhaps with large length scales and a 
large magnetic helicity might lead to an equilibrium. The range of 
equilibria in an isentropic star is however more restricted than that 
available in a stably-stratified star, and the possibility that an equi- 
librium could be reached from any realistic initial conditions seems 
rather unlikely. In addition, the convective motions driven by heat 
flux, which are not included in these simulations, would presum- 
ably prevent equilibrium formation. These questions will be studied 
more fully in a forthcoming publication. 



3.3 Ongoing flux loss during accretion 

We have seen how flux can be destroyed once the protostar has 
stopped accreting and is fully detached from its surroundings, but 
it seems plausible that something similar to these two processes re- 
move flux in accreted material as, or soon after, it is accreted. This 
is possible because unless the field is already extremely weak the 
Alfven timescale on which these mechanisms work is much shorter 
than the accretion timescale. However, since we do not observe the 
star until the main accretion phase is over, it is difficult to distin- 
guish observationally between ongoing flux destruction during the 
embedded accretion phase and flux destruction only at the end of 
that phase when the star becomes detached. 

In the ma gnetospheric accretion model dKoenigll 1 19911; 
lLong et al]|2008[) . there is a gap (with low plasma /3) between the 
star and the inner edge of the disc; gas is channelled from the inner 
edge along field lines anchored around the magnetic poles of the 
star. The bulk of the stellar flux self-connects within the magneto- 
sphere, so the stellar field is essentially detached from the surround- 
ings - it is implicitly assumed in this model that some mechanism 
exists to break field lines between the star and the parent cloud. 

In addition, note that the field strength drops off sufficiently 
quickly with increasing distance from the star that the magnetic 
energy in the vicinity of the star is much greater than that further 
away. [More precisely, the magnetic energy per unit radius as a 
function of radius (£ , ma g('") = B 2 r 2 /2) must drop off faster than 
log normal, meaning that dlnB/dlnr < —3/2; the values for 
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Figure 6. A time sequence at t/r\ = 0, 1.6, 4.2, 9.4, 13.1 and 19.3 (top row left to right and then bottom row left to right). The field quickly reaches the 
surface, but only in some places, since there must be a downwards flow in other places. Discrete flux tubes become recognisable, as in the fourth frame. In 
the last frame, the entire field can be thought of as discrete flux tubes. As these flux tubes are pushed upwards by buoyancy they lose their axial component 
into the atmosphere which causes them to become longer and thinner. This is similar to process in stably stratified stars, except that it happens on a dynamic 
timescale rather than a diffusive timescale. 



dipole and split-monopole fields are —3 and —2 respectively.] This 
means that energy liberated by reorganisation of the field inside 
and immediately outside the star outweighs any energy required to 
change the geometry of the field further away. 

In the general case where a net non-zero (but probably small) 
helicity has been accreted, a radiative star must contain an equilib- 
rium which evolves quasistatically as a result of the addition of new 
mass, flux and helicity. Imagine that the star has a weak dipole field 
and accretes material with magnetic field aligned to that dipole: the 
material can simply flip over until its field is antiparallel to that of 
the star (similarly to the case of two bar magnets in section UTt 
whereupon its field is annihilated against that of the star; the field 
lines originally connecting the blob to the outside world will re- 
connect so as to bypass the star. In this way a lower energy state is 
reached, magnetic energy being converted into heat. In fact, since 
the flux of the accreted material is not conserved, only the accreted 
mass and helicity are relevant for the evolution of the star and its 
magnetic field. 

If on the other hand the star is convective or has a significant 
convective envelope, its field will be in a dynamo-powered steady 
state. Highly magnetised material arriving at the surface will stay 
at the surface until both its entropy and magnetic buoyancy can be 



reduced by mixing and reconnection. Again, the accreted flux is 
not relevant, but here also the accreted helicity is not necessarily 
conserved, depending on the properties of the dynamo. 

If these mechanisms can prevent significant flux from building 
up in the star, the accretion disc will not be able to drag magnetic 
flux inwards; the magnetic field lines necessarily find some way 
of slipping outwards relative to the inspiralling gas. This may be 
facilitated by an outwards-directed Lorentz force resulting from the 
build up of flux towards the centre. When the disc first forms it must 
initially be able to drag flux inwards, accounting for the large flux 
observed in the inner disc; at some point that inwards drag switches 
off or is much reduced as magnetic pressure builds up in the centre. 



4 DISCUSSION 

Let us assume for the moment that a significant fraction of the 
original flux survives transport into the protostar. In low-mass 
(< OAMq) protostars which have not yet become fully convective, 
it is possible that the flux is destroyed while accretion continues, 
but at the very latest when it ceases. If it is possible to measure the 
magnetic field on a low-mass pre-MS star after significant accretion 
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has ended but before deuterium ignition, one should see a global 
MHD equilibrium rather like those seen in radiative stars on the 
MS. It is probably difficult or impossible to measure the magnetic 
field on low-mass protostars which are still accreting significantly, 
but this measurement could distinguish between ongoing flux de- 
struction and destruction only after complete detachment from the 
surroundings. In principle flux destruction could begin as soon as 
the Alfven timescale becomes less than the collapse timescale, i.e. 
when the first core forms and the collapse becomes sub-Alfvenic. 
Once deuterium ignites however, and the star becomes fully con- 
vective - which it remains throughout the main-sequence - the star 
is expected to lose all of its original flux and helicity and with it the 
'magnetic initial conditions' from the parent cloud. This explains 
why the magnetic fields of low-mass stars appear to depend only on 
rotation speed. Solar-type stars (eventual mass 0.4 — 1.5Mq) also 
pass through a fully-convective phase before a radiative core devel- 
ops; in this envelope a dynamo similar to that in low-mass stars is 
expected, and observed. Inside the radiative core, one might expect 
to find the remnant of the previous convective dynamo in an equi- 
librium whose strength depends on the helicity generated by (or 
surviving) this dynamo. Note that for a dynamo to generate a net 
helicity from zero initial helicity there must be symmetry-breaking 
or positive feedback; alternatively an initial net helicity inherited 
from the molecular cloud may persist through the dynamo phase if 
the phase does not last too long, or may perhaps be amplified. In 
any case, the residual field is unlikely to be stronger than equipar- 
tition with the previous convective energy. 

Intermediate-mass stars (> 1.5M©) also have a fully- 
convective protostellar phase , but convection in the envelope even- 
tually dies away. In the radiative envelope we are therefore seeing 
either formerly convective material, which has lost all 'memory' of 
its initial magnetic field, as in the solar-type stars, or material which 
accreted onto the star after significant convection disappeared, in 
which case the accreted helicity may survive. Thus the difference 
between magnetic field strengths in intermediate-mass MS stars 
and the bimodality between magnetic and non-magnetic stars could 
be at least partially due to differences in the accretion history; the 
non-magnetic stars would have accreted all of their material while 
still convective, while the magnetic stars would have become ra- 
diative while accretion continued and were able from that point on- 
wards to accumulate magnetic helicity from the accreted gas. This 
is possibly connected to the observational result that the magnetic 
fraction i ncreases with mass, from ~ 1% at 1.5M@ to ~ 20% 
at 6Mq jPowej 2007). In more massive stars deuterium fusion 
contributes relatively little to the energy and convection may be 
unimportant, in which case we do expect to see magnetic equilibria 
which depend on initial conditions. However, if a star forms from a 
relatively symmetrical cloud then helicity should be small and we 
expect almost all of the magnetic energy to be lost during the relax- 
ation to equilibrium. An equilibrium can then survive essentially 
unchanged for the entire main-sequence, owing to the long Ohmic 
timescale (> 10 10 yr), although it is possible that rotation-induced 
circulation has some effect over main-sequence timescales. 

An important unsolved question therefore is on what timescale 
can a convective dynamo generate or destroy helicity. It is plausible 
that a dynamo approaches a true steady state (erasing initial condi- 
tions) on the diffusive timescale associated with the length scale 
of the convection (Brandenburg, priv. comm.) which may be rather 
longer than a pre-MS convective phase, but probably shorter than 
the ages of low-mass MS stars. 

Another unsolved problem is that of the enormous range 
in field strengths seen in stars with radiative envelopes. It has 



been suggested (Zinnecker & Yorke 20071: 


Maitzen et al.1 120081; 


Boeomazov & Tutukov 2009; Ferrario et al. 


20091) that magnetic 



stars are merger products, the merger resulting in differential rota- 
tion and some kind of dynamo activity. This scena rio would explain 
the la ck of magnetic stars in short-period binaries dAbt & Snowdenl 
Il973h . However, this could also be explained by magnetic inhi- 
bition of fragmentation, as is seen in simulations dMachida et al.l 
2008); also there is the issue of the lack of universal magnetic fields 
in blue stragglers. In any case, it seems likely that the magnetic stars 
are unusual in some sense and that the normal state of affairs is for 
the flux to drop to that corresponding to less than 1 gauss. 



5 CONCLUSIONS 

I have examined mechanisms to destroy flux in a protostar, with 
a view to shedding some light on the observational fact that the 
magnetic flux in main-sequence stars is very small compared to 
that in molecular cloud cores in star-forming regions. It is uncer- 
tain how much of the original flux is already removed from the 
gas by various diffusive processes during collapse and accretion, 
and the mechanisms explored here can remove flux which survives 
these processes and is accreted onto the protostar. The mechanisms 
can be categorised according to whether they operate in convective 
or non-convective protostars. In convective stars, the magnetic field 
rises buoyantly to the surface of the star on a dynamical timescale 
(an Alfven timescale, ~ 10 yr in a solar-type star with a 1 kG field 
or somewhat longer in a protostar with larger radius) and its en- 
ergy is destroyed by reconnection in the atmosphere, much as is 
observed in the solar corona, until the field strength has dropped to 
that which can be maintained by a convective dynamo. In a radia- 
tive star (or the radiative zone of a star), an arbitrary magnetic field 
evolves on the same dynamical timescale into an equilibrium, the 
strength of which depends not on the original flux but on the mag- 
netic helicity. A non-zero helicity requires some asymmetry in the 
magnetic field of the accreted material; in the standard hourglass 
model the helicity is zero. 

In summary any excess 'unwanted' flux can be destroyed once 
the star becomes magnetically independent from its surroundings, 
via buoyancy and/or MHD instability and reconnection on a dy- 
namical timescale. This should happen either while the main ac- 
cretion phase is ongoing or as it comes to an end - in any case, by 
the time we can observe the star it will already have lost its original 
flux. 
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